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We extend the application of vector and axial Ward identities to calculate 6a, bp and br, coefficients that give 
the mass dependence of the renormalization constants of the corresponding bilinear operators in the quenched 
theory. The extension relies on using operators with non-degenerate quark masses. It allows a complete deter- 
mination of the O(a) improvement coefficients for bilinears in the quenched approximation using Ward Identities 
alone. Only the scale dependent normalization constants Z% (or Zg) and Zt are undetermined. We present 
results of a pilot numerical study using hadronic correlators. 



To remove errors of 0(a) from physical matrix 
elements, one must improve both the action and 
the operators jy. The former requires the addi- 
tion of the Sheikholeslami-Wohlert term 101 , 



Ssw = -a 5 c S w '^'^{x)j<J lll ,F lll ,{x)il}{x) . (1) 



Improvement of flavor off-diagonal bilinears re- 
quires H both the addition of extra operators, 

{M)n = A^ + ac A d^P (2) 
(Vi) M = V^ + acycUV (3) 
(Tj)^ = T M!/ + acr(d„V v - d^) , (4) 

and the introduction of a mass dependence, 

= z£(l + & x am y )(Xj)W>. (5) 

Here X = A,V,P,S,T, Z\ are the renormal- 
ization constants in the chiral limit, and rriij = 
(rrii + rrij)/2 are the bare quark masses defined 
using the axial Ward Identity (WI), Eq. (g). The 

bare unimproved bilinears are Atf = ^7^75 ijj^ , 
etc. The task is to determine the coefficients csw, 
Z^'s, ex's, and &x's non-perturbatively. 

Previous calculations have shown how Z v , Z\ 
and Z°/Z s §, c SW , c A and b v §§, c v g, 
ct ||, and bp — b A and 65 ^1 can ^e deter- 
mined non-perturbatively using axial and vector 
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WI. We discuss here an extension that yields b A , 
bp, and bx- The two remaining constants Z P (or 
Zg) and Zip are scale and scheme dependent, and 
so cannot be determined using WI. Note that the 
relations we derive do not extend directly to the 
unquenched theory, which requires additional im- 
provement constants and a more complicated set 
of conditions . 

We begin by recalling the ALPHA method for 
determining csw and c A |(|. The bare WI mass 

9 _ Eg(W + ac A d^Pp\x,t)J^\0)) 

should be independent of t, up to corrections 
of 0(a 2 ), since it is proportional to the average 
rcnormalized quark mass 



to, + mf 
1 j 



Z A (1 + b A arriij) 
Z p (l + bparriij) 



(7) 



This is achieved by simultaneously tuning csw 
and c A . Our approach differs from the 
Schrodinger functional method of Ref. || in that 
we use standard 2-point correlation functions. 
Consistency of these estimates are checked by 
varying the initial state using J = P or Ai and 
with different types of sources for the quark prop- 
agators (Wuppertal smearing, Wall, point). In 
the following we use the abbreviation rrii = ma, 
where ma refers to two degenerate flavors. 



2 



With c$w fixed, Z v and by are obtained using 
charge conservation. We use the forward matrix 
elements of (Vi)4 between pseudoscalars, 



Z v (l + byam 2 ) 

(^P(i2)(f,r)J(2D(0)) 



(8) 



with r > t > and J = P or A 4 . Note that the 
cy term in Vj does not contribute. 
Next consider the generic axial WI 



SS^0f\y)J^(z)) = (S0^(y) jW(z)) (9) 
where O 23 = ^ 2 )l\A (3) , SO 13 = t/jW^TipW, and 



12) 



(mf + m«)(P fi )( 12 )-a Al (A i? )( 12 )l(10) 



This results from a chiral rotation on flavors 1, 2 
in the 4-volume V, with y € V and z ^ V. By 
enforcing these identities in the chiral limit one 
can determine cy and cx HJ8|^| j as shown below. 

Away from the chiral limit, operators P and O 
in the product / v (mf + m«)(P i? )( 12 )(a;)0^ 3) (y) 
need off-shell improvement. This requires the ad- 
dition of a contact term, of unknown normaliza- 
tion, having the same form as the RHS of (^|) ||] . 
Our new observation is that the contact term is 
proportional to mi+m2 and so can be removed by 
extrapolating mi and m,2 to zero. This leaves the 
freedom to examine the dependence on 7713, and 
from this one can determine certain combinations 
of the bx- In the following, the extrapolation to 
mi = t«2 = is implicit^]. 

As a first application of this method we show 
how to obtain 6a, as well as cy, using the AWI 



' - i 



Z° A {l + b A am 3 /2) 



Z\ ■ Z v (l + b v am 3 /2) 

E ? <(^/)i 13) (£s/4) J< 31 )(0)) 

E e^(6SW ( Vl )f 3) (y, y4 )4 (31) (0)) 



E ? e^((^)f 3) (^y4)Ar^(o)) 



(31), 



(11) 



(12) 



2 In practice, we keep the f (mi + »i2)(-Pfl)' 12 ' term in Eq. 
( |lo[ ) prior to extrapolation, since it improves the signal. 



where J = P or A4,. Eq. (|ll| ) is independent of 
cy &tp= 0, and its m 3 dependence gives bA — by. 
The intercept provides a second determination of 
Z v . Eq. ( |l2] ) is used to determine cy. 

Given cy, an alternate determination of bA — by 
is obtained from 



r-2 



Z v (l + b v am 3 /2) 
Z° A ■ Z%(1 + b A am 3 /2) 
E e ^-^ 5 (12) {Al )f 3 \y,y 4 ) 



V, 



(13) 

(3D (0)) 



y/ 31) (0)) 



This also yields Z A . The same information can 
be obtained from the combinations 



1 



- 7° 

— A A 1 



r 2 

Z° 

- = ^{l + {bA-by)am 3 /2) 



(14) 
(15) 



[n 

Similarly, we determine bp — bg and Z p /Z s using 
Z P {1 + b P am 3 /2) 



7 n 



Z s (l + b s am 3 /2) 

J2 ff e^y(P^)(y 7 y 4 ) J( 3 D(0)) ' 
To get c T we use the WI with O = 



(16) 



1 + acT 



.(31) 



(31) 
k4 



(0)> 



(17) 



At p = 0, {Ti)ij has no contribution from the Ct 
term, so the only cx dependence is on the LHS. 

The previous method fails for bx since both 
sides of (||) have the same dependence on bx if 
mi = m 2 — 0. The cure is to consider three 
non-degenerate masses. The contact term re- 
quired to improve the LHS of (||) is proportional 
to mi + ni2, while the relative dependence of the 
two sides on bx is proportional to mi — ra-i- Thus 
the two terms can be separated and bx deter- 
mined. More details will be given in |11|. We 
have not yet implemented this proposal. 

Thus far we do not have a separate determina- 
tion of bp or bg. This can be accomplished with 
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the method of Ref. 10 1, which uses 2-point cor- 
relation functions for non-degenerate masses. In 
addition to using that method we present a vari- 
ant which avoids the need to study quantities as 
a function of the underlying hopping parameter 
k. We first note that if we use Eq. (0) and enforce 
(2mi) R + (2m 2 ) R = 2(mi + m 2 ) R , we find 

4m i2 - 2 [mn + TO22] 



b A 



(18) 



a[m n - m 22 \« 
We next make use of the vector two-point WI 



vmi2 



^ e v-(,S( 12 )(f,t)J( 21 )(0)) 
Z v [l + b v a(m 1 +m 2 )/2 



Z%[1 



Am i2 



-t bsa(mi + m 2 )/2] 

where the source J is either </( 21 ) = S^ 21 * 1 or 
^_. P (23)(- Z4 )p(3i)( ) fo r < t < 24- Enforc- 
ing 2(mf — m^) = (2mi) fl — (27712)^, we find 
bs - b v 



(bp -I 
Ami2 



R z [m 
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m 22 \ 



aRz[m[ 1 — m\ 2 \ 

where i? z = ZP S Z\I{Z%Z%). Since 6^ and 6 y are 
already known, 6p and bs are given by Eqs. (18) 
and ©. 

The above discussion shows that in principle 
one can determine all the constants, except Z P 
(or Zg) and Z^., in the quenched theory using 
Ward identities. The results of an exploratory 
study are summarized in Table jj] These were ob- 
tained on 83 lattices of size 16 x 48 at (3 = 6.0. 
Since the action is only tree-level tadpole im- 
proved (csw — 1-4755), the results do not rep- 
resent full 0(a) improvement, but they indicate 
the efficacy of the method. More details will ap- 
pear in 1 1 1 1 . 

We draw two preliminary conclusions. First, 
even though we have found channels in which the 
statistical and systematic errors on the determi- 
nation of b A — by and bp — bs are fairly small, the 
magnitude of these differences are still compara- 
ble to their error. This rough equality between all 
the bx is consistent with perturbative results jL2| . 
Second, the determination of cy has a large un- 
certainty, which accounts for a substantial frac- 
tion of the errors in Z A , Z P /Zg, and ct- For ca 



Eq.# 


observable 


intercept 








C A 


-0.016(11) 




g 




z% 


+0.746(1) 








V 


+ 1.55(2) 


( 


11 


) 


zi 


+0.752(7) 


( 


11 


) 


bA — by 


+0.34(21) 


( 


12 


) 


cy 


+0.46(29) 


( 


14 


) 


a VI \ A) 


+1 32(12") 


( 


14 


) 


bA — by 


+ 1.8(1.1) 


( 


14 


) 


Z A 


+0.78(2) 


( 


15 


) 


zyz v 


+ 1.00(5) 


( 


15 


) 


bA - by 


+ 1.2(8) 


( 


16 

IT) 


) 


^A^Sl^P 


+0.96(1) 


( 


) 


bp - b s 


-0.08(9) 


( 


17 


) 


ct 


-0.14(7) 


1 


10 
10 


1 


70 7U 17V) 
^A^sl^P 


+0.96(1) 


I 


1 


b A -b P + b s /2 


+0.49(1) 


( 


18 


) 


bp - b A 


+0.1(4) 


( 


19 


) 


bs-by- 2(b P - b A ) 


-0.5(5) 



Table 1 

Constants extracted from the different WI. 

and cy, the Schrodinger functional method j| || 
gives results with much smaller errors, and may 
prove to be the method of choice. 
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